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ABSTRACT 


The  effects  of  sound-speed  and  current  variations  induced  by  a 
mesoscale  cyclonic  eddy  on  short-range  propagation  are 
considered.  A parametric  eddy  model  is  used  to  determine 
acoustically-relevant  eddy  environmental  effects,  so  that  eddy- 
acoustical  effects  can  be  determined  for  eddies  of  arbitrary 
size,  strength,  and  position.  Approximations  to  sound-speed 
and  current  structures  are  used  to  investigate  eddy  effects 
on  the  three-dimensionality  of  rays  and  on  ray  types.  The 
influence  of  current  and  sound-speed  variations  on  travel  time 
is  examined,  and  accurate  expressions  for  per-ray  phase  variation 
are  obtained.  Examples  are  presented  illustrating  effects  of 
source-receiver  position  and  orientation  on  per-ray  phase  shifts 
and  relative  phase  spreading  of  arrivals.  Also,  general  results 
are  presented  which  illustrate  the  variations  of  eddy-acoustical 
effects  as  functions  of  source-receiver  range  and  of  eddy  size 
and  strength. 
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INTRODUCTION 

Recent  studies  of  the  mesoscale  structure  of  some  mid-ocean 

regions  have  shown  that  eddies  are  often  common-place  and  long- 

lasting  features.1  These  large  rotating  bodies  of  water  have 

been  shown  to  have  significant  effects  on  the  sound-speed 

2 

structure  of  the  ocean,  primarily  resulting  from  the  large 

temperature  deviations  associated  with  eddies.  Several  analytical 

studies  of  the  environmental  influences  of  eddies  have  been 

carried  out,  including  Refs.  3 and  4.  The  authors  of  this  paper 

have  presented  a model  of  acoustical ly-relevant  sound-speed 

perturbations  and  currents,  relating  these  environmental  effects 

to  parameters  characterizing  eddy  size  and  strength.'* 

Some  studies  have  demonstrated  that  eddies  can  have  substantial 

effects  on  the  transmission  of  sound  in  the  ocean.  Nearly  all 
2 6 7 

investigators  ' ' have  assumed  two-dimensionality  of  acoustic 

rays  in  studying  propagation  through  a sound-speed  structure 

computed  numerically  from  oceanographic  measurements  of  an  eddy. 

8 9 

However,  Baer  ' has  considered  both  two-  and  three-dimensional 
parabolic-equation  methods  to  investigate  propagation  through  an 
eddy,  using  our  model  in  Ref.  5.  Each  of  the  above  studies, 

though,  considered  only  a particular  eddy  of  fixed  size  and  strength, 
and  no  investigation  was  made  of  the  dependence  on  eddy- 

acoustical  effects  of  eddy  variations.  Further,  these  studies 
did  not  consider  arbitrary  source  and  receiver  orientation  with 
respect  to  the  eddy  involved.  Moreover,  they  ignored  the  effects 
of  the  large  rotational  currents  found  in  eddies,  which  have  been 
suggested  to  have  possibly  significant  acoustical  effects. ^ 
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The  objective  of  this  paper  is  to  present  general  results 
for  the  influence  of  ocean  eddies  of  arbitrary  size  and  strength 
on  the  transmission  of  sound.  Both  currents  and  sound  speed  are 
included,  and  no  restriction  is  placed  on  the  horizontal  location 
of  the  acoustic  source.  In  this  initial  analytical  study, 
however,  a short  range  (<50  km)  problem,  with  both  source  and 
receiver  on  the  surface,  is  considered.  Future  studies  by  the 
authors  will  permit  larger  ranges,  and  other  source  and  receiver 
depths.  In  this  paper,  an  investigation  qf  both  two-  and  three- 
dimensional  effects  on  ray  geometry  is  carried  out,  and  eddy 
effects  on  per-ray  travel  time  and  phase  are  determined  as 
functions  of  eddy  parameters  and  source  and  receiver  position. 

In  Sec.  I,  approximations  to  both  eddy  currents  and  sound- 
speed  effects,  for  use  in  the  acoustical  transmission  problem, 
are  derived.  Using  a previously-derived  parametric  eddy  model,  these 
are  expressed  as  simple  functions  of  eddy  size  and  strength. 

Section  II  then  investigates  the  effect  of  eddy  sound-speed 
variations  and  currents  on  ray  geometry,  employing  an  environmental 
discretization  scheme.  The  three-dimensionality  of  rays  is 
examined,  and  found  to  be  a higher-order  effect  in  short-range 
problems.  The  types  and  variation  of  rays  within  an  eddy  are 
studied,  and  ray  paths  determined  for  a fixed  source  and  receiver. 
Section  III  investigates  eddy  effects  on  a transmitted  acoustic 
signal.  The  principal  effects  are  found  to  arise  through  per-ray 
travel  time,  and  hence  phase  variation,  and  analytical  expressions 
are  obtained  for  these  quantities.  Several  examples,  relating 
source  and  receiver  position  and  orientation  within  an  eddy  to 
per-ray  phase  variation,  are  considered  in  Sec.  IV.  Also,  current 
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effects  are  shown  to  be  significant,  so  that  they  should  not  be 
ignored  in  eddy-acoustic  studies.  General  results,  illustrating 
the  utility  of  an  analytical  procedure,  are  presented  for  the 
maximum  variation  of  phase  as  a function  of  range  and  eddy 
parameters.  It  is  shown  that  phase  variations  of  over  20  cycles 
are  possible  for  each  ray.  Section  V summarizes  the  principal 
results  of  the  paper. 

I.  SOUND-SPEED  AND  CURRENT  MODELS 

As  stated  previously,  the  authors  have  developed  one  model 

for  sound-speed  and  current  variations  induced  by  a class  of 

mesoscale  eddies.^  After  appropriate  non-dimensionalization 

and  scaling  of  the  fluid-dynamical  equations,  approximate 

solutions  for  currents  and  motion-induced  pressure  and  density 

perturbations  were  obtained.  Resulting  temperature  effects  were  then 

determined  from  the  Eckart  state  equation, ^ after  which  sound  speed 

1 2 

was  found  from  a modified  Wilson's  formula.  These  quantities 

were  taken  to  be  radially  symmetric  about  the  eddy  center,  with 
horizontally-uniform  vertical  decay  and  with  dependence  on  eddy 
radius,  maximum  current  speed,  and  depth  of  influence.  The  two 
coordinates  of  the  eddy  center  completed  the  parameter  set  for  our 
model . 

To  determine  sound-speed  effects  of  an  eddy  of  given  size 
and  strength,  perturbation  density  and  pressure  must  first  be 
computed  from  the  model  equations  and  used  in  conjunction  with 
prescribed  static  density  and  pressure  in  the  Eckart  and  Wilson 
formulas.  Typical  resulting  sound-speed  profiles  at  the  eddy 
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edge  and  center  are  shown  as  the  solid  curves  in  Fig.  1(a) 
for  a large  cyclonic  eddy  of  radius  125  km  and  maximum  current 
100  cm/sec.  The  principal  effects  are  a decrease  in  sound 
speed c within  the  eddy,  resulting  from  the  cold  core,  and  a 
significant  elevation  of  the  SOFAP.  axis  of  over  300  m.  As  a 
result  of  the  complexity  of  the  Wilson  and  Eckart  formulas,  the 
dependence  of  these  effects  on  both  location  within  the  eddy  and 
model  parameters  is  nonlinear.  However,  we  shall  utilize  bilinear 
approximations  to  the  sound-speed  structure  to  simplify  both 
this  functional  dependence  and  ray  geometry. 

Let  D denote  the  uniform  depth  of  the  ocean,  and  z the 

s 

SOFAR-axis  (minimum  sound-speed)  depth  in  the  static  case  - that 
with  the  eddy  absent.  We  take  the  surface  and  bottom  sound 
speeds  equal  and  denote  them  by  cq,  which  will  be  taken  to  remain 
constant  with  or  without  the  eddy  present.  The  minimum  sound 
speed  at  the  SOFAR  axis  in  the  static  case  is  represented  by 
c . Specification  of  D,  z.  , c , and  c completely  determine 
the  bilinear  approximation  to  the  static  sound  speed,  as 
illustrated  by  the  lower  dashed  curve  in  Fig.  1(a).  In  order 
to  model  the  eddy-induced  variations,  we  first  specify  a 
Cartesian  coordinate  system,  centered  on  the  surface  at  the  eddy 
center.  The  x and  y coordinates  lie  on  the  ocean  surface,  while 
z is  measured  positively  downward.  The  SOFAR  axis  depth  at  any 
horizontal  position  (x,y)  is 

ZA  - zAs  11  - <la) 

where 


the  SOFAR-axis  depth,  and  6 (r)  = 0 outside  the  cylindrical  eddy 


(r  > r ) . We  note  from  Fig.  1(a)  that,  below  the  SOFAR  axis 


sound-speed  variations  are  minimal,  so  that  the  lower  sound-speed 


constant  value 


Variations  in  sound  speed  above  the  axis  are  represented  by  a 


so  that  the  sound-speed  minimum  c is  given  by 


and  the  corresponding  upper-gradient  magnitude  is 


We  note  that  the  perturbed  sound-speed  structure  at  any  point 


within  the  eddy  is  then  determined  through  Eqs.  (1)  - (3)  by 


Again,  due  to  the  nonlinearity  of  the  Eckart  and  Wilson 


formulas,  it  is  not  possible  to  derive  a simple  analytical  form 


However,  numerical  calculations  have  been  made  for 
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static  parameters  appropriate  to  the  mid-ocean  Atlantic  region. 
These  have  shown  that  the  maximum  SOFAR-axis  elevation,  ^zAmax' 
can  be  written  as  a simple  function  of  the  eddy  radius  rQ  and 
maximum  current  speed  UQ,  and  is  relatively  insensitive  to  the 
depth  of  influence  of  the  eddy.  From  Fig.  9 of  Ref.  5,  it  is 
seen  that  this  relationship  is  almost  linear  in  both  Uq  and  rQ, 
so  that  we  may  write 


Aza  = 0.028  U r , 
Amax  o o 


(4) 


where  U is  measured  in  cm/sec,  r in  km,  and  Az„  in  m. 
o o Amax 

Further,  the  model  equations  show  that  the  sound-speed  deviations, 
and  hence  the  fractional  SOFAR-axis  elevation  6(r),  have  the 
following  properties: 


and 


6(0)  = max  6 (r) , 


Mro)  = 0 , 


6'  (0)  =6*  (rQ)  = 0 . 


(5a) 

(5b) 

(5c) 


We  select  the  following  form  for  6 (r)  that  is  motivated  by  the 
model  equations  and  matches  the  solution  exactly  for  small  axis 

elevat ions : 


L 


Mr)  = (flrpjnait/rAsHJ0(or/ro)  - JQ  (a)  ][1-J0  («)  ] 


-1 


(6) 


In  Eq.  (6),  JQ  is  the  Bessel  function  of  the  first  kind  of  order 
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zero,  and  a = 3.83  is  the  first  zero  of  J^,  the  Bessel  function 
of  the  first  kind  of  first  order.  Equations  ( 1 ) — ( 6 ) completely 
specify  the  sound-speed  approximations  at  any  point  within  an 
eddy  of  known  radius  and  maximum  current  speed. 

Similar  approximations  are  now  made  to  the  current  structure 
within  an  eddy.  The  maximum  current  speed  at  any  horizontal 
location  is  at  the  surface,  z = 0,  and  decays  monotonically  to  zero 
at  the  maximum  depth  of  eddy  influence,  z=z  , as  illustrated  by  the 
solid  curve  in  Fig.  1(b)  for  a surface  current  magnitude  of 
100  cm/sec.  Separability  of  eddy  vertical  and  horizontal 
structure  in  the  model  equations  implies  that  the  rate  of  this 
decay  is  identical  at  any  horizontal  location,  so  that  only  the 
magnitude  Sq  of  the  surface  current  varies,  given  from  the  model 
equations  as 


SQ(r)  =(UQ/m)  J L (ar/rQ) , 


where  m =0.582  is  the  maximum  value  of  J^.  We  then  approximate 
the  vertical  current  decay  from  its  surface  maximum  by  a linear 
decrease  to  zero  at  the  SOFAR  axis,  while  neglecting  the  small 
currents  below  the  axis.  The  resulting  expressions  for  the 
current  components  u and  v in  the  x and  y directions , respectively , 
are  then  linear  functions  of  depth  (see  Fig.  1(b)),  have  surface 
magnitude  SQ,  and  depend  on  position  in  the  eddy  as  follows; 


( u , v ) = sQ(r) [ (zA~z)/zA) (-y/r,  x/r) . 
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II.  RAY  GEOMETRY  AND  ENVIRONMENTAL  DISCRETIZATION 

We  consider  the  case  of  a sound  source  and  a receiver  & 

both  on  the  ocean  surface  a distance  R apart.  For  convenience, 
we  orient  our  coordinate  system  with  origin  at  the  eddy  center  and 
with  x-axis  parallel  to  the  source-receiver  direction.  Using 
this  geometry,  the  source  is  located  at  (xs,yg) , a distance  rg 
from  the  eddy  center,  and  the  receiver  at  (xR,  yR)  with  yR  = yg, 
a distance  r from  the  center,  as  illustrated  in  Fig.  2(a). 

x\ 

As  a result  of  both  eddy  currents  and  horizontal  sound-speed 

variations,  an  acoustic  ray  launched  from  the  source  and 

intercepting  the  receiver  will  not  remain  in  the  vertical  plane 

containing  d and  (%  ; rather, it  will  be  a three-dimensional 

space  curve.  To  examine  the  deviation  out  of  a vertical  plane,  we 

consider  a ray  launched  initially  in  this  plane  toward  the  receiver. 

. 1 3 

One  of  the  modified  scalar  ray  equations  proposed  by  Ugincius 

A 

describes  the  deviation  y from  the  xz-plane.  After  discarding 
lower  order  effects  arising  from  vertical  currents  and  products 
of  small  quantities,  this  equation  simplifies  to 

(cqc  1 y')  = z ' cq  ^ 3v/3z  - co*9c/3y  - x'cq^  9u/3y  , (9) 

where  primes  denote  differentiation  with  respect  to  arc  length, 
and  (x',  y',  z')  are  the  direction  cosines  of  the  rays.  The 
terms  on  the  right  side  of  Eq.  (9)  describe  the  effects  that  can 
cause  the  rays  to  become  three-dimensional.  The  first  results 
from  the  vertical  gradient  of  the  current  component  v perpendicular 


- 


V 


k 
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to  the  source-receiver  direction.  The  second  is  associated 
with  the  horizontal  sound-speed  gradient  out  of  the  source- 
receiver  plane,  while  the  third  involves  the  variation  of  the 
current  component  u in  the  source-receiver  direction.  Estimates 
of  the  maximum  magnitudes  of  these  gradients,  obtained  numerically 
from  figures  in  Ref.  5,  are  |3v/3z|  = 0(10  sec  ), 

|3c/3y|  = 0(10  ^sec  ^)  , and  |3u/3y|  = 0(10  ~*sec  ^).  Thus  current 
effects,  represented  by  the  first  term  on  the  right  of  Eq.(9), 
can  be  more  influential  in  causing  three-dimensional  effects  than 

9 

previously  studied  sound-speed  effects. 

Equation  (9)  should  be  integrated  twice,  along  with  two 
additional  equations  describing  the  x-  and  z-  behavior  of  the 
rays.  However,  to  the  order  of  terms  retained  in  Eq . (9) , it 

can  be  shown  that  it  is  valid  to  perform  these  integrations  along 
the  static  two-dimensional  paths  connecting  source  and  receiver. 

A 

To  obtain  an  upper  bound  for  the  excursion  y of  a ray,  forms  for 
u,  v,  and  c were  obtained  from  Eqs.  (1)  - (8) , with  the  use  of 
large,  but  reasonable,  parameter  values.  For  ranges  less  than 
50  km,  which  are  the  only  ones  to  be  considered  in  this  paper,  the 

/v 

deviation  y was  found  to  be  less  than  50  m.  This  leads  to  a 
horizontal  bearing  error  of  less  than  0.01  deg,  and  an  error 
of  less  than  0.1  m in  total  arc  length,  resulting  in 
insignificant  errors  in  phase.  Hence,  we  will  assume  that  all 
rays  are  two-dimensional.  For  longer  ranges,  however,  deviations 
out  of  the  xz-plane  can  be  of  the  order  of  hundreds  of  meters, 
leading  to  significant  bearing  errors  and  phase  variations,  so 
that  the  three-dimensionality  of  rays  cannot  necessarily  be 


10 


neglected.  Results  for  larger  source-receiver  separations, 

utilizing  horizontal  sound-speed  variations,  have  been  obtained 
9 

by  Baer.  However,  only  one  particular  eddy  was  considered, 
and  eddy  currents  were  not  included. 

A typical  idealized  ray  from  source  to  receiver  is  shown 
in  Fig.  2(b).  We  index  this  ray  by  the  number  N of  bottom  reflections 
and/or  bottom  turning  points  strictly  between  and  </-^  , so 
( that  the  figure  corresponds  to  an  N = 3 ray.  We  note  that,  for 

l 

the  short  range  (R  < 50  km)  problems  considered  here,  the 
distance  between  bottom  reflection  and/or  turning  points  is  small 
compared  to  the  eddy  radius.  Thus,  little  variation  in  eddy 
structure  is  observed  over  this  distance.  Therefore,  we  approximate 
the  eddy  sound-speed  and  current  structure  by  taking  average 
vertical  profiles  between  the  source  and  the  first  bottom  point, 
between  bottom  points,  and  between  the  last  bottom  point  and  the 
receiver.  This  partitions  the  range  R into  N+l  intervals, 

as  illustrated  by  the  dashed  lines  in  Fig.  2(b),  and  in  each 
interval,  current  and  sound  speed  are  functions  of 
depth  z only.  From  Eqs.  (1)  - (8)  this  discretization  can  be 
accomplished  by  defining  the  mean  values  of  the  fractional  SOFAR- 
axis  elevation,  6^,  and  the  surface  current  speed  in  the  x direction. 
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- • - ~ ■ -r  .. 


X . 

i = -<X1+1  - X,,-1  J1 


[SQ(r)ys/r]dx 


(10b) 


We  note  that  the  integration  is  taken  on  the  surface  of  the 

ocean,  along  part  of  the  line  connecting  the  source  and  receiver, 

and  that  the  values  of  for  each  ray  are  unknown  at  this  time. 

We  now  consider  a ray  launched  in  the  xz-plane  from  the 

source  at  an  angle  0 measured  positive  clockwise  from  the 

horizontal . As  a result  of  the  horizontal  discretization  in 

Eqs.  (10),  in  each  interval  both  the  current  u in  the  source- 

receiver  direction  and  the  sound  speed  c are  functions  of  depth 

alone,  given  by  Eqs.  (2),  (3),  (8),  and  (10).  For  cases  of 

purely  depth-dependent  sound-speed  and  current  structures, 

14 

Franchi  and  Jacobson  have  shown  that  the  angle  0 which  the 
ray  makes  with  the  horizontal  at  any  point  can  be  related  to  the 
initial  launch  angle  0qo  through  a modified  Snell's  law: 


-1  a -2  -1  _ -2 
c cos  0 - c V = c cos  0 - c u. 

o oo  o oo 


The  angle  0 . at  the  surface  in  interval  I.  is  given  by 
oi  1 


cos  0 . = cos  0 + (V  . - V ) c , 

Ol  OO  Ol  oo  o 


(12a) 


while  the  angle  0 = 0^  at  the  bottom  point  x^  is  given  by 


cos  0,  . = cos  0 - V c 

bi  oo  oo  o 


(12b) 


I 
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If  the  right  side  of  Eq.  (12a)  is  greater  than  one  in  any  interval, 
the  ray  will  not  reach  the  surface,  and  it  becomes  what  is  commonly 
called  a "mixed"  ray.  Similarly,  if  the  right  side  of  Eq.  (12b) 
is  greater  than  one,  the  ray  will  refract  before  reaching  the 
bottom  and  also  become  mixed.  However,  since  the  current  terms 
relative  to  cq  in  Eqs.  (12)  are  small,  a ray  will  only  become 
mixed  when  the  initial  angle  0oo  is  small.  The  largest  possible 
current  magnitude  in  the  eddy  is  the  surface  maximum  U • If  0QO 
is  sufficiently  large  so  that 

cos  0 + U c ^ < 1 , (13a) 

oo  o o ' 

the  ray  will  remain  SRBR  (surf ace-ref lected/bottom- ref lected) 

at  all  points  within  the  eddy.  For  small  0 and  U c \ this 
c 1 oo  o o 

inequality  can  be  approximated  by 

0 > (180/it)  (U  c-1)1/2  , (13b) 

oo  o o 

where  Qqq  is  measured  here  in  degrees.  For  a maximum  eddy  current 
speed  Uq  of  150  cm/sec,  all  rays  launched  at  angles  0qo  > 3.6° 
remain  SRBR.  However,  numerical  calculations  show  that,  for 
ranges  less  than  50  km,  all  rays  that  intersect  the  receiver  have 
a launch  angle  greater  than  6°  a>nd,  hence,  remain  SRBR.  We 
stress  that,  as  range  is  increased,  current-induced  mixed  rays 
cannot  be  neglected,  and  their  existence  will  greatly  influence 
the  received  signal. 

14 

It  has  been  shown  that,  when  sound-speed  and  currents  are 
only  linear  functions  of  depth,  ray  paths  can  be  approximated  by 


» 
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whore  zA.  is  Eq.  (la)  evaluated  at  6 = 6..  Similarly,  c . is 
Ai  l J mi 

the  axial  sound-speed  from  Eq.  (3a)  when  <5  = <5^.  In  a manner 

analogous  to  Ref.  15,  it  can  be  shown  that  the  distance  d to 

n 


the  n th  surface  reflection,  for  a ray  launched  at  the  angle 


14 


The  angles  0 ^ and  0^  at  the  surface  and  bottom,  respectively, 
are  given  in  Eqs.  (12) , while  the  angles  0m^  at  the  SOFAR  axis 


cos  0 . = c . c cos  0.  . 
mi  mi  o bi 


from  Eq.  (11).  Examples  of  d.  and  the  angles  0 0,_  . , and  0 . 

1 3 oif  bi'  mi 

appear  in  Fig.  2(b). 

To  determine  the  launch  angle  of  a ray  that  intersects  the 
receiver  with  N bottom  reflection  points,  Eq.  (15a)  is  set  equal 
to  the  range  R.  The  expressions  for  R^ f ! ®mi'  ®Qi'  Rci^  ' an<^  ®bi 
from  Eqs.  (12),  (14),  and  (16)  are  substituted,  and  the  result 

is  simplified  further  by  using  expressions  for  g. . , z .,  V . and 
c . . It  is  then  necessary  to  determine  the  launch  angle  0 both 
to  completely  specify  the  ray  geometry  via  Eq.  (11)  and  to  permit 
determination  of  variation  of  quantities  such  as  spreading  loss 
and  bottom  and  surface  phase  and  amplitude  effects.  However,  the 
complexity  of  the  resulting  equation  precludes  exact  solution  for 
the  launch  angle  0 ._  It  is  possible  to  solve  for  this  angle 

in  the  static  case,  0QS»  where  no  currents  or  horizontal  sound- 
speed  variations  are  present: 


tan  0 = (c  + c ) (DN)  (c  R)  1 - R(c  -c  ) (4Nc  D)  1.  (17a) 

os  o ms  o o ms  o 


We  then  seek  an  approximation  to  0qo  by  writing 


tan  0 = tan  0 + y. 

oo  os 


(17b) 
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We  note  that  Eq.  (17b)  is  appropriate  when  the  quantities  in  Eq.  (15a) 
differ  only  slightly  from  their  static  counterparts.  The 
approximation  implied  by  Eq.  (17b)  is  thus  accurate  near  the  eddy 
edge  where  the  SOFAR-axis  elevation  remains  small.  However,  near 
the  eddy  center,  elevations  in  axial  depth  occur  that  are  large  enough 
to  produce  6 = 0(10 _^),  so  that  significant  errors  are  conceivable 
in  using  Eq . (17b). 

To  avoid  this  problem,  we  define  a new  environmental  state, 

where  variations  along  the  transmission  range  remain  small. 

This  "barred"  state  has  no  currents  or  horizontal  sound-speed 

variations.  Moreover,  the  sound-speed  structure  at  all  points 

is  identical  to  that  at  the  source  location,  with  axial  depth 

z.  and  sound-speed  minimum  c , given  by 
A m 


ZA  ■ zAS[1-5lrS>' 


(18a) 


and 


g2  zAs  MrS>' 


(18b) 


The  barred  sound-speed  state  is  an  exact  representation  at  the 
source,  and  if  range  is  small  compared  to  eddy  radius,  sound-speed 
variations  over  the  transmission  range  are  small.  The  amounts 
of  the  variations  can  again  be  defined  through  a fractional  axis 
elevation  <5^  in  interval  1^.  In  terms  of  the  6^,  we  can  express 
the  axial  depth  and  sound-speed  minimum  at  any  point  as 


z 


Ai 


Z A (1'  V 


(19a) 


\ 


and 
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c . = c - g_  z.  6.  , 
mi  m 2 A l 


(19b) 


where  is  given  by 


5.  = (6.  - 6 (rs)] [1  - 6 (rs)] 


(19c) 


from  Eqs . (la)  and  (18a). 

As  deviations  from  the  barred  state  are  small  for  transmission 
ranges  of  interest  here,  we  seek  the  launch  angle  6qo  in  the  form 


tan  8 = tan  8 + y , 

oo  o 


(20a) 


where  8q  is  the  launch  angle  in  the  barred  state,  given  analogously 
to  Eq.  (17 a)  by 


tan  6 = (c  + c ) (DN) (c  R)  1 - R(c  -c  ) (4c  ND)  1 . (20b) 

o o m o o m o 


Equation  (20a)  can  be  used  in  Eq.  (15a)  (with  dN  = R)  to  determine 
the  launch  angle.  Using  expressions  for  currents  from  Eqs.  (7)  and 
(8) , as  well  as  the  sound-speed  deviations  from  the  barred  state 
implied  by  Eqs.  (19),  the  resulting  equation  may  be  simplified. 

It  will  then  depend  only  upon  ocean  parameters  and  deviations 
from  the  barred  state,  represented  by  cq^  , 6^,  and  y . These 

last  quantities  are  now  sufficiently  small  that  squares  and 
products  can  be  neglected,  further  simplifying  the  equation  for 
8oq.  After  lengthy  algebraic  simplification,  the  perturbation  to 
the  launch  angle  is  determined  in  the  form 


• ...  f t ■ * 
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v = *!  £ + »2  V^c;1  + a3  £ M(i,voico1  ' (21a) 


where 


zA[cos  (8m~0o) -1] [ 2N (d-za) cos  0Q(sin  0m  - sin  0Q)  ] • (21b) 


2-  - -1 
- [cos  0 sin  0 ] 
o o 


(21c) 


z.(l  - sin  0 sin  0 ) sin  0 [2ND  cos^8  sin  0 (sin  0 -sin0  )]  ^ 
A o m m oo  mo 


- zR  [2ND(cos  0 - cos  0 )] 

A o m 


In  Eqs.  (21),  the  quantity  cos  8^  is  obtained  from  Snell's  Law: 


cos  6 = c c cos  0 . 

m m o o 


(21d) 


(21e) 


The  first  term  on  the  right  in  Eq.  (21a)  results  from  sound-speed 
variations  from  source  to  receiver,  the  second  from  the  current 
at  the  source,  and  the  third  from  variations  in  the  current  from 
source  to  receiver.  We  note  that  a^,  a 2,  and  a^  are  all  0(1), 
so  that  the  correction  y to  the  launch  angle  is  of  the  order  of 
the  small  parameters  v0^c~^  and  6^. 

The  magnitude  of  the  angle  perturbation  y depends  on  the  sum 


of  the  6.  and  V ..  From  Eqs.  (10), 
l oi 


E Mfi)*i = E M(i)  !(xi+i_xirl  ( 

• /»  • ** 


x . . - 
i*i  + l 


6 ( r ) dx  - 6 (rs) j [1-6 (rg) ] 


x . . 
,1  + 1 


X>U)Voi  = £ M(i)|(xltl-xi)-1  ( S0(r)y 

• * f\ 


.r  ^dx!  . 


(22a) 


(22b) 
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The  bottom  reflection  points  are  still  undetermined,  since 
the  launch  angle  0qq  must  be  known  to  completely  describe  the 
geometry  of  a ray.  However,  by  examining  an  equation  similar  to 
Eq.  (15a)  for  the  location  of  these  points,  it  can  be  shown  that 
if  the  sound-speed,  current,  and  angle  terms  vary  only  by  terms  of 
0(y),  then  the  fractional  offset  (x^  - x^)/x^  is  0(y).  Hence, 
consistent  with  neglecting  products  of  small  parameters  in  the 
derivation  of  Eq.  (21a),  x^  + ^ - x^  can  be  replaced  by  x^  + ^ - 
The  relationship 


xi+1  - xi  = M(i) R(2N) 


(23) 


can  then  be  used  to  simplify  Eqs. 
N 

X)  M(i)6.  = 2N[l-<5  (rg)  ]_1  Jr1 
i=0 


to 


6 (r) dx 


and 


6(rS) 


(24a) 


N 

E 


2NR 


-1 


S" 


So(r,*S 


-1 


dx 


(24b) 


1 = 0 XS 

where  the  integration  is  taken  along  the  surface  of  the  ocean 
on  a line  connecting  source  and  receiver. 

The  correction  y to  the  static  angle,  necessary  to  determine 
the  variation  in  boundary-interaction  effects  and  spreading  loss, 
can  be  expressed  from  Eqs.  (17b)  and  (20a)  as 


y = tan  0 - tan  0 + y . 

' o os 


(25a) 
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Further  simplification  is  possible  using  Eqs.  (17a)  and  (20b) , 
resulting  in 

Y = Y + Kco  - cms)  ZAsl  ( C0  ( D - zAs)l 
x [ DNR~ 1 - R(4ND)_1]6 (rg)  , 

where  y is  given  in  Eq.  (21a).  Numerical  computations  show  that 

for  the  ranges  of  under  50  km  considered  here,  the  launch-angle 

variation  0 - 0 is  less  than  one  degree.  Hence,  from  Eq.  (11), 

oo  os 

the  corresponding  angle  variation  at  any  other  point  also  satisfies 
this  bound.  However,  for  longer  ranges  where  the  launch  angle 
becomes  smaller,  the  presence  of  sine  terms  in  the  denominators  of 
Eqs.  (21)  could  cause  much  larger  angle  variations. 

The  approximate  ray  geometry  is  now  completely  specified  for 
a fixed  source  and  receiver.  As  a result  of  the  horizontal  sound- 
speed  and  current  discretizations,  all  rays  are  approximated  by 


segments  of  circular  arcs.  In  summary,  for  short-range  problems, 
we  have  shown  that  all  rays  may  be  taken  to  be  two-dimensional, 

. 

and  that  only  SRBR  rays  intercept  the  surfaced  receiver.  These 
rays  are  indexed  by  the  number  N of  bottom  reflection  points, 

| ' I 

beginning  with  the  N=1  ray.  Through  use  of  the  small  variation 
in  eddy  structure  over  range,  an  accurate  expression  for  the 
launch  angle  from  Eq.  (20a)  can  be  obtained  for  each  ray.  This 
is  accomplished  by  calculating  0Q  and  y from  Eqs.  (20b),  (21), 

and  (24).  This  procedure  is  valid  for  arbitrary  eddy  size  and 
strength,  as  well  as  for  arbitrary  source-receiver  orientation 
and  position  within  the  eddy. 


i 


u Mwmmm 


- 20  - 

III.  TRAVEL  TIME  AND  PHASE  VARIATION 

Specification  of  all  rays  connecting  source  and  receiver  and 
their  geometries  permit  determination  of  travel  times.  From 
Ref.  14,  travel  time  T may  be  computed  as 

T = u cos  0)  1 ds,  (26) 

where  the  integral  is  taken  along  the  ray  path  from  to  . 

However,  Hamilton,  Siegmann,  and  Jacobson^  have  shown  that 

contributions  to  the  travel  time  due  to  deviations  in  ray  paths, 

caused  by  environmental  variations  of  the  order  of  some  small 

2 

parameter  e,  are  of  higher  order,  0(e  ),  than  explicit  contributions 
from  the  integrand.  Such  an  argument  has  been  used^  to  then 
evaluate  Eq.  (26)  along  the  known,  simpler  static  path.  However, 
in  this  study,  environmental  deviations  from  the  static  ocean 
may  be  large,  of  0(6),  so  that  significant  errors  might  arise  by 
neglecting  squares  of  these  terms.  This  difficulty  may  again  be 
avoided  by  exploiting  the  "barred"  state.  As  the  barred  sound-speed 
state  is  exact  at  the  source,  environmental  variations  such  as 
currents  and  sound-speed  variations  will  be  small  (0(y)] 

Consistent  with  our  previous  neglect  of  products  of  these  small 
parameters,  ray  path  variations  can  be  ignored  in  evaluating 
Eq.  (26).  Consequently,  the  integral  in  Eq.  (26)  may  be  taken 
along  the  corresponding  ray  path  that  connects  the  source  and 
receiver  in  the  barred  state,  where  sound  speed  is  only  a function 
of  depth  and  no  currents  are  present. 

Further  approximations  can  also  be  made  to  the  integrand 
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in  Eq.  (26) . It  is  consistent  to  expand  the  integrand  in  the 
small  quantities  representing  sound-speed  variations  and 
currents  scaled  by  the  surface  sound  speed  cq,  and  to  neglect 
their  powers  and  products.  Thereby,  Eq.  (26)  can  be  simplified  to 

T = |{[c(z)J  1 - Ac(x,z)c  2 - u(x,z)c  2 cos  0 } ds,  (27a) 

J o o o 

where  the  sound  speed  c has  been  written  as 

c = c(z)  + Ac(x,z).  (27b) 

In  Eqs. (27)  c(z)  is  the  sound  speed  at  the  source,  the  "barred"  state 
given  in  Eqs.  (2),  (3),  (6),  and  (18),  and  Ac(x,z)  is  a convenient 
representation  of  sound-speed  variation  from  source  to  receiver  as 
described  in  Sec.  II. 

To  determine  eddy  effects  on  travel  time  and  phase,  it  is 
advantageous  to  write  T as  the  static  travel  time  Tg  plus  small 
eddy-induced  fluctuations  AT: 

T = T + AT.  (28) 

s 

From  Eqs.  (26)  and  (27),  AT  can  be  expressed  as 

AT  = (T  - T ) + AT  + AT  , (29a) 

s c u ' 

where  T is  the  travel  time  in  the  barred  state. 


22 


E 


r i 

■.i 

It: 


\ 


■ 


t 


at 

c 


/ 


- I Ac  (x, z) c 


-2 


ds 


(29b) 


is  the  contribution  from  sound-speed  variations  between  source  and 
receiver,  and 


at 

u 


/ 


u (x, z) c 


-2 

o 


cos  0 ds 
o 


(29c) 


results  from  eddy  currents.  From  Ref.  15  the  travel  times  in 
the  static  and  barred  states  can  be  evaluated  as 

(30a) 

(30b) 


T = [ 2ND (c  -c  ) log[c  c } (1+sin  0 ) (1+sin  0 ) ] 

s o ms  ^ o ms  ms  os 

and 


T = [2ND(c  -c  ) log[c  c 1 (1  + sin  0 ) (1  + sin  0 ) 1 ] , 

o m o m m o 


where  log  denotes  the  natural  logarithm-  The  difference  (T  - Tg) 

in  Eq.  (29a) , representing  the  variation  in  travel  time  between 

the  true  static  case  outside  the  eddy  and  the  barred  statef  can  be 

simplified  first.  Because  the  magnitudes  of  T and  Tg  given  in 

Eqs.  (30)  are  very  similar,  it  is  most  accurate  to  simplify 

their  difference  by  first  combining  the  logarithm  terms.  The 

resulting  expression  can  then  be  linearized  in  parameters 

characterizing  differences  in  sound  speed  between  the  two  states 

-2 

divided  by  cq.  These  parameters  remain  small  [0(10  )], 

although  the  difference  in  axial  depth,  5 ( ) , can  be  large  [O(10_1)]. 
Before  writing  the  final  expression  for  T - Ts  , we  turn  to  the 
remaining  terms  in  Eq.  (29a). 
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To  evaluate  the  integrals  in  Eqs.  (29b)  and  (29c),  the  forms 
for  the  sound-speed  variations  Ac  and  currents  u from  Sec.  II  are 
used.  The  ray  paths  are  connected  circular  arcs,  so  that 


ds  = 


c (g  cos  0 ) 1 d0 
o 1 o 


z < z. 


co(g2  COS  V 


-1 


de 


z > z 


Ola) 


and 


z = 


z. (cos  0 - cos  0) (cos  0 -cos  0 ) ^ 

A o o m 

z +(D-z  ) (cos  0-cos  0 ) (cos  0 -cos  0 ) ^ 

A A mom 


A ' 


7 > 7 

z - A 

z > z.  . 
A 


(31b) 


From  Eqs.  (3),  (10)  and  (19),  the  sound-speed  variations  Ac.  and 


currents  u^(z)  in  range  interval  1^  are  given  approximately  by 

( -g.^  D 6^  (D-za)  1 2 0 < z < zA 

0 


Ac . = 

l 


(32a) 


z > z 


A 


and 


V .(z-z)z-  ^ 
oi  A A 


u . = 

l 


0 < z < z. 


z > z . 


(32b) 


Again,  products  of  terms  that  will  arise  in  Eqs.  (29b)  and  (29c) , 
such  as  fractional  SOFAR  elevation  and  fractional  sound-speed 
variation,  are  neglected. 

After  extensive  manipulation  of  the  expressions  for  T - T^ 
and  the  integrals  in  Eqs.  (29b)  and  (29c)  (using  Eqs.  (31)  and 
(32)],  the  travel-time  change  AT  may  be  written  as 

XR  *R 

AT  = Ksl  6(rS)+Ks2[1_<S(rS)  1 J Vq^tJc”1  dx+Ks3  J [6  (r)-6  (rs)  ]dx.  (33) 


In  Eq.  (33),  K . , K _,  and  K depend  only  on  the  static  state 

S -L.  S A S j 

of  the  ocean,  and  may  be  written  as 


„ , T _ . r m ( 1 - cos(0  - 0 ) 

Ksl  = 2nDzAs  Co  D_ZAs  > 08  ms 


(cos0  -cos9  ) (sin0  -sin0  ) 
os  ms  ms  os 


(a^d+a^)  1cos0qs]  + a2[a1(l+a1)  ^-log  ( 1+a.^ ) +log  ( l-a2 ) +a2  ( l-a2 ) * ] j , 

(34a) 


s2 


, ( (0  -0  )cos0  - (sin0  -sin0  ) ) 

, „ r\  . v-1  ) ms  os  os  ms  os  ( 

'As  1 ' os  co  j (cos0  - cos0  ) (sin@  - sin0  )(  ' 

\ os  ms  ms  os  > 


and 


K = [z (c  -c  ) ]( (D-z )c2]  1 

s3  As  o ms  As  o 


sin0  -sin0  -(0  -0  )cos0  \ 

ms os ms  os ms  / 

(cos0  -cos0  ) (sin0  -sin0  )( 
os  ms  ms  os  I 


(34c) 


In  Eqs.  (34) 


and 


a,  = R(cos  0 ) (c  -c  ) [2ND(l+sin  0 )c  ] 

1 os  o ms  os  o 


-1 


-1 


= (c  -c  )c 
2 o ms  o 


<34d) 


( 34e) 


We  note  that  all  of  the  quantities  in  Eqs.  (34)  are  determined 
independent  of  the  size  and  strength  of  the  eddy  and  the  positions 
of  the  source  and  receiver  in  the  eddy.  That  is,  the  coefficients 
Ksr  Kg2>  and  Kg^  in  the  per-ray  travel-time  variation  are  determined 
entirely  by  the  transmission  range  R and  the  number  of  bottom 
reflections  N [since  0 and  0 are  determined  by  these  quantities 
through  Eqs.  (11)  and  (17a)].  The  size  and  strength  of  the  eddy. 


- 25  - 


as  well  as  the  orientation  of  the  source  and  receiver,  enter 
in  the  determination  of  AT  in  Eq.  (33)  solely  through 
the  axial  variations  6(r)  and  the  surface  currents  VQ(r). 

The  variation  AT  requires  computation  of  integrals  of 
SOFAR  axis  elevations  6 (r)  and  surface  current  components,  VQ(r), 
from  the  source  to  receiver.  We  note  that  closed-form  evaluation 
of  these  integrals  is  impossible,  due  to  the  complex  form  of 

6 and  VQ  given  in  Eqs,  (6)  — (8)  . However,  we  choose  to  exploit  the 
fact  that  our  range  R is  small  compared  to  the  eddy  radius  rQ , so 
that  the  integrals  are  approximated  by  the  trapezoidal  rule: 


A _ 

f 

) [6 (r)  - 6<r 


s)]dx  = R [ 6 ( rR)  - 6(rs)]/2 


(35a) 


f - 

J V (r)  c 


o (r) c0  dx  = R(VQ(rR)  + Vq ( rg) ] / ( 2cq) . 


(35b) 


Our  range  assumption  also  can  be  used  to  simplify  the  forms  of 

K . , K 9 , and  K ...  Since  all  rays  are  relatively  steep,  the 
S X S « S 

. . —2 
quantity  0ms  - 0os  is  small  [0(10  )].  By  expanding  the 

trigonometric  terms  in  Eqs.  (34b)  and  (34c)  in  powers  of  ®ms“®os 

and  discarding  higher  order  terms,  it  can  be  shown  that 


zAs(2Dco) 


z.  [c  -c  ][2(D-z.  )c2cos0  ] ^ 
As  o ms  As  o os 


(36a) 


(36b) 


26 


f I 

I i 
t i 


These  expressions  are  most  accurate  for  ranges  shorter  than 
approximately  30  km,  but  still  represent  lowest-order  effects 
for  larger  ranges.  The  expression  for  can  also  be  simplified 

by  considering  propagation  between  a source  and  receiver  both 
lying  on  the  same  diameter,  with  the  source  on  the  eddy  edge. 

As  all  currents  will  be  perpendicular  to  an  eddy  diameter,  no 
current  effects  enter  into  the  determination  of  AT  through  Eq.  (33). 
It  follows  from  Eqs.  (33)  and  (35a)  that  the  path  reciprocity  of 
travel-time  variations  if  source  and  receiver  positions  are 
reversed  implies 


si 


= R K 


s 3 ‘ 


(37) 


We  note  that  a corresponding  expression  exists  if  the  trapezoidal 
rule  approximation  is  not  applied,  and  that  Eq . (37)  has  been 

approximately  verified  numerically.  Equations  (35) -(37)  can  then 
be  used  to  simplify  Eq . (33)  giving 


AT  = ! Rz.  (c  -c  )[2(D-z.  )c2  cos0  ] 
I As  o ms  As  o os 


-1 


- R Z-  ( 2D  c2)  1 [l  - 6 ( r ) ] [V  (r  ) 


o ' S' 


(6(rs)  + 6 (rR) ]/2 
* Vo(rR),/2  ' (38) 


The  first  term  on  the  right  of  Eq.  (38),  corresponding  to  the 
sound-speed  variation  within  an  eddy,  represents  the  effect  of 
colder  water  found  in  a cyclonic  eddy,  in  which  a travel  time 
delay  will  occur.  This  delay  depends  only  on  the  the  cosine  of 
the  initial  launch  angle  (or,  equivalently,  the  transmission  range) 
for  any  particular  ray.  Since  the  cosine  decreases  with  the 
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number  of  bottom  reflections,  larger  phase  shifts  will  occur  for 
each  successive  arrival.  Thus,  ray  arrival  times  will  spread 
as  a result.  The  second  term  in  Eq.  (38)  gives  the  effects  of 
eddy  currents  on  the  phase  of  each  ray.  For  a current  component  in 
the  same  (opposite)  direction  of  propagation,  a ray  will  travel 
faster  (slower) , and  a negative  (positive)  travel-time  shift  will 
occur.  We  note  that  there  is  no  dependence  in  this  term  on  the 
number  of  bottom  reflections,  so  that  the  current-induced  shift  is 
uniform  over  all  rays.  This  result  represents  a balance  between  two 
effects;  specifically,  the  additional  time  that  rays  with  more 
bottom  reflections  spend  in  the  eddy  and  are  exposed  to  eddy 
currents,  and  the  decrease  with  N of  the  current  component  in  the 
ray  direction,  u cos  0q. 

The  preceding  determination  of  travel-time  deviations  for 
each  ray  permits  ready  construction  of  per-ray  phase  shift.  We 
consider  an  omnidirectional  CW  source,  emitting  a signal  with 
frequency  f(in  Hz).  The  phase  <f>N  (in  cycles)  of  a ray  having  N 
bottom  reflections  may  be  written  as^ 

<}>N  = fT  + (N-l)/2  + NSn  . (39) 

The  first  term  on  the  right  represents  the  contribution  of 
travel  time;  the  second,  that  of  a 7T-rad  phase  shift  at  each  of 
the  N-l  surface  reflections;  and  the  third,  that  of  the  phase  shifts 
(in  cycles) at  each  of  the  N bottom  reflections.  We  note  that  eddy 
effects  are  assumed  to  enter  Eq.  (39)  only  through  the  travel 
time.  This  is  because  the  numbers  of  surface  and  bottom  reflections 
remain  unchanged  here,  and  the  bottom  phase  shift  SN  is  relatively 


insensitive  to  small  changes  [0(y)l  in  incident  angle,  for  the 


steep  rays  characteristic  of  a short-range  problem.  Therefore 


where  ATM  is  the  travel -time  change  from  Eq.  (38) 


The  amplitude  of  each  ray  A^  can  be  determined  by  finding 
spreading  loss  and  bottom  loss  per  reflection.  However,  the 


variation  of  the  ray  angles  at  each  bottom  reflection  and  at 


the  receiver  have  been  shown  to  be  small  [0(y)l.  It  follows,  as 


order  importance  in  determining  the  ray  amplitude.  Hence,  the 


easily-computed  static  spreading  loss  and  bottom  loss  expressions 


If  desired,  the  amplitude  A and  phase  0 of  the  total  acoustic 


field  at  the  receiver  may  be  found  from  the  sum  of  all  ray  arrivals 
at  , each  of  which  has  an  amplitude  A^^  and  phase  <t>N.  However, 
in  mid-ocean  regions,  clay-silt  bottom  conditions  most  often 


which  have  been  shown  to  lead  to  high  bottom  losses 


for  the  steep  angles  characteristic  of  short  range  problems 


Under  such 


bottom  conditions 


decreases  rapidly  with  the  number  of  bottom  reflections,  so  that 


the  N=1  ray  constitutes  the  principal  contribution  to  the  total 
acoustic  field.  It  is  then  possible  to  approximate  the  total- 
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little  variation  in  A will  occur.  The  total-field  phase  <!>  can 
also  be  approximated  accurately  by  the  phase  <t>^  of  the  N=1 
ray,  so  that  variations  in  total-field  phase  can  be  determined 
from  Eqs.  (39)  and  (38)  with  N=l.  VJe  do  stress,  however,  that 
if  a less  lossy  bottom  is  assumed,  it  is  necessary  to  consider  the 
effects  of  additional  rays.  In  such  a case,  the  phase  spreading 
of  consecutive  arrivals  may  lead  to  large  total-field  amplitude  and 
phase  fluctuations,  often  characteristic  of  multipath  effects. 

IV.  APPLICATIONS 

We  now  discuss  application  of  the  results  of  our  preceding 
theory  to  some  particular  types  of  sound  transmission  problems. 

We  stress  the  generality  of  these  results.  This  generality  will 
be  a consequence  of  our  analytic  treatments  of  eddy  current  and 
sound-speed  effects  for  arbitrary  eddy  radius  r , maximum  current 
speed  UQ,  range  R,  and  source  and  receiver  position  and 
orientation,  given  by  (xg,  yg)  and  (xR,  yR) . From 
Eqs.  (6)  - (8),  (17a),  (38),  and  (40),  the  explicit  dependence  of 

the  change  of  phase  (in  cycles)  for  the  N th  ray,  due  to  the 

presence  of  an  eddy,  is 

A<|>  = f]  B U r (R2  + B^N2  + B^R4N_2  + 2R2B  B ) 1/2 

N (loo  2 3 23 

x (JQ(ars/ro)  + JQ(arR/ro)  - 2JQ(a)] 

+ B4  R U0[J1(ars/ro)',SrS1  + Jl(t,rR/ro)',RrRl1 

* (1  ' B5Vo'V“rS/ro>  ' V°>l>!  ’ <4J«> 

where 


In  Eq.  (41a),  the  dimensional  quantities  B^ , B2,  B^ 
only  on  static  ocean  parameters  and  are  given  as 

B1  - °-°28  <co-cmS>  !2  co(D'zas> 


(41b) 


(41c) 

and  B^  depend 


(42a) 


= (c+c) D 
o ms 


(42b) 


(42c) 

( 4 2d) 
(42e) 


In  Eqs . (41)  and  (42),  currents,  sound  speeds,  depths,  and  radial 

distances  are  measured  in  units  of  cm/sec,  m/sec,  m,  and  km, 
respectively.  Equations  (41a)  can  then  be  used  to  illustrate 
variation  of  the  phase  of  each  ray,  as  a function  of  the  parameters 
characterizing  eddy  size  and  strength  and  source-receiver  range, 
position,  and  orientation. 

It  is  first  necessary  to  specify  the  quantities  in  Eqs.  (42) 
by  specifying  the  static  state  of  the  ocean  area  under  consideration. 
Here,  we  elect  to  use  representative  data  taken  from  the  mid-Atlantic 
ocean,  as  illustrated  by  the  solid  curve  in  Fig.  1(a) . ^ A 
uniform  depth  of  5 km,  surface  sound  speed  cQ  = 1543  m/sec, 
and  a static  SOFAR-axis  depth  z^g  = 1100  m,  with  a corresponding 
sound-speed  minimum  cmg  = 1483  m/sec,  reasonably  fit  the 
static  sound-speed  profilers  illustrated  by  the  lower  dashed 
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curve  in  Fig.  1(a).  We  also  consider  a source  emitting  a CW 
frequency  f =400  Hz.  Specification  of  these  values  gives 
numerical  values  of  = 6.46  * 10  = 9.82  * 10^, 

B3  = 1.94  * 10  6,  B4  = 7.94  x lcT8,  and  B5  = 1.81*  1CT5,  These 

parameter  values  will  be  used  for  all  calculations  in  this  section. 

We  first  investigate  the  effect  of  variation  of  source 
and  receiver  position  within  an  eddy.  The  eddy  radius  rQ  is 
fixed  at  125  km  and  the  maximum  current  speed  UQ  at  100  cm/sec, 
values  typical  of  a Gulf  Stream  Ring.^  The  maximum  SOFAR- 
axis  elevation  is  obtained  from  Eq . (4)  to  be  350  m.  Initially, 

the  source  and  receiver  are  taken  to  be  outside  the  eddy,  on 
the  ocean  surface,  separated  by  a range  R = 30  km.  The  source- 
receiver  line  is  then  moved  through  the  eddy  center,  perpendicular 
to  an  eddy  diameter,  as  illustrated  in  Fig.  3.  Our  choice  of 
coordinate  system  implies  that  this  diameter  may  be  taken  in  the 
direction  of  the  y-axis  with  the  receiver  moving  on  the  diameter. 
Such  a path  could  also  result  from  an  eddy  drifting 

appropriately  across  a fixed  transmission  range,  We  note 
that,  in  following  such  a path,  the  receiver  will  move  through 
the  point  of  maximum  sound-speed  variation  at  the  eddy  center. 

It  will  also  pass  through  points  of  maximum  surface  currents, 
at  about  half  a radius  above  and  below  the  center. 

Since  eddy  radius,  maximum  current  speed,  and  transmission 
range  are  fixed,  the  only  variation  of  the  phase  shift  of 
each  ray  in  Eqs.  (41)  enters  via  source  and  receiver  position. 

At  each  point  along  the  source-receiver  path,  source  and  receiver 
position  can  be  determined,  and  the  per-ray  phase  shifts 
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calculated  through  Eq.  (41a) . These  are  illustrated  in  Fig.  3 
for  the  N=l,2,  and  3 rays.  Initially,  both  the  source  and 
receiver  are  outside  the  eddy,  resulting  in  the  static  phase, 
and  hence  a phase  variation  of  zero.  As  they  move  into  the  eddy, 
the  slower  sound  speeds  in  the  cold  core  result  in  phase  increases 
which  reach  a maximum  at  the  eddy  center.  Rays  having  more 
bottom  reflections  spend  more  time  in  the  eddy,  and  hence  experience 
a larger  phase  increase,  resulting  in  the  relative  spreading  of 
arrivals  shown  by  the  vertical  distances  between  curves  in  Fig.  3. 
Effects  of  eddy  currents  are  illustrated  by  the  dashed  lines  in 
the  figure.  For  negative  (positive)  y,  below  (above)  the  eddy 
center,  currents  are  in  (opposed  to)  the  direction  of  propagation, 
so  that  a smaller  (larger)  phase  increase  occurs.  Hence,  each 
phase  curve  is  not  symmetric  about  y = 0,  as  it  would  be  if  currents 
were  neglected.  The  horizontal  dashed  lines  indicate  a maximum 
difference  of  about  one  cycle  in  phase  shift  observed  at  equal 
distances  on  either  side  of  the  eddy  center. 

A similar  behavior  would  be  observed  if  the  source  and 
receiver  were  to  follow  any  other  chord,  not  necessarily  a 
diameter,  through  an  eddy.  Sound-speed  effects  would  cause  the 
phase  shifts  of  each  ray  to  increase  as  the  source  and  receiver 
move  toward  the  center  of  the  eddy,  and  then  decrease  after 
passing  the  point  of  closest  approach  to  the  center.  Currents 
would  increase  the  amount  of  phase  shift  when  opposite  to  the 
direction  of  propagation,  and  decrease  the  phase  shift  when 
flowing  from  source  to  receiver.  We  note  that  current  effects 
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are  represented  by  the  second  term  on  the  right  side  of  Eq.  (41a) . 

As  discussed  earlier,  the  current  effects  are  independent  of  Nf 
so  that  this  variation  is  uniform  for  all  rays.  However,  the 
magnitudes  of  current  and  sound-speed  effects  will  vary  with  the 
path  of  the  chord.  For  example,  if  the  transmission  range  fails 
to  intersect  the  eddy  center,  the  maximum  observed  phase  shifts 
will  be  smaller. 

The  orientation  of  the  source  and  receiver  at  any  point  will 
also  influence  the  observed  phase  shifts,  particularly  with 
respect  to  current  effects.  We  again  fix  eddy  radius  at  125  km  and 
maximum  current  speed  at  100  cm/sec.  A receiver  is  fixed  60  km 
above  the  eddy  center,  and  the  source  is  rotated  clockwise  about 
the  receiver,  keeping  a fixed  range  of  30  km,  as  illustrated  in 
Fig.  4.  Thus, only  the  orientation  and  radial  distance  of  the  source 
vary.  Again,  the  only  information  required  to  calculate  eddy- 
induced  phase  shift  in  Eq.  (41)  is  the  position  of  the 
source  and  receiver.  The  resulting  phase  variation  of  the  N=1 
ray  is  plotted  in  Fig.  4 as  a function  of  rotation  angle  6 
relative  to  the  arbitrary  initial  orientation  described  by  6=0. 

The  dashed  curve  illustrates  the  effect  of  sound-speed  variations 
alone.  That  is,  currents  are  absent,  as  in  the  case  of  previous 
eddy-acoustics  studies.  The  behavior  of  the  curve  illustrates 
the  initial  motion  of  the  source  into  an  area  of  lower  sound 
speed,  towards  the  eddy  center,  resulting  in  larger 
(positive)  phase  shifts.  For  6 > 180°,  the  source  has  moved 
into  an  area  of  higher  sound  speed,  toward  the  eddy 
edge,  resulting  in  smaller  (negative)  phase  shifts.  The  larger 
variations  at  90°  as  compared  to  270°  occur  because  higher 
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sound-speed  variations  arise  moving  toward  the  eddy  center  than 
toward  the  eddy  edge. 

The  solid  curve  in  Fig.  4 illustrates  the  significant  phase 
variations  produced  by  inclusion  of  currents.  At  the  arbitrary 
initial  angle  0=0,  currents  are  entirely  in  the  direction  of 
propagation.  As  the  source  is  rotated,  current  components  in 
the  source-receiver  direction  first  decrease  and  then  become 
negative,  until  at  0 = 180°  the  current  is  opposite  to  the 
direction  of  propagation.  Initially,  the  current  component  at 
the  source  increases  while  that  at  the  receiver  decreases,  so 
little  variation  is  noted  as  0 increases  to  about  30°.  Further 
increases  in  0 result  in  larger  phase  variations  than  would 
result  if,  as  in  previous  studies,  only  sound-speed  effects 
were  included.  The  largest  variation  occurs  at  0 = 180°,  where 
an  error  of  one  cycle  arises  if  currents  are  neglected,  as 
shown  by  the  vertical  arrow  in  Fig.  4.  This  is  compared  to  a 
total  variation  of  5 cycles  if  only  sound-speed  effects  are 
considered.  Again,  as  indicated  by  Eq.  (41a),  current  effects  will 
be  identical  for  other  rays.  At  any  other  location  in  an  eddy, 
results  similar  to  those  in  Fig.  4 would  be  observed.  A variation 
in  sound-speed  effects  would  occur  as  the  source  moves  closer  toward 
or  further  from  the  eddy  center.  Similar  current-induced  variations 
would  occur  as  the  current  component  in  the  source-receiver 
direction  changes  with  rotation  angle.  The  magnitudes  of  these 
effects  would  vary  with  the  magnitudes  of  sound-speed  variations 
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and  currents,  as  indicated  in  Eq.  (41a).  Both  of  the  examples 
illustrated  in  Figs.  3 and  4 clearly  indicate  the  necessity  of 
considering  the  effects  of  eddy  currents  on  the  phase  of  an 
underwater  signal. 

A principal  advantage  of  an  analytical  treatment,  as  opposed 
to  a numerical  study,  is  that  general  results  can  be  obtained. 
Thus,  through  Eqs.  (41),  the  effects  on  acoustic  phase 
variations  of  transmission  range,  eddy  radius,  and  maximum 
current  speed  can  be  examined  easily.  We  choose  to  discuss 
effects  of  these  variations  on  the  maximum  per-ray  phase 
variation,  which  results  from  the  large  sound-speed  variations 
near  the  eddy  center.  Indeed,  by  examining  Eqs.  (41)  for  small 
source  and  receiver  radii,  it  can  be  shown  that  the4 largest  per- 
ray  phase  variation  (denoted  by  A<J>N)  occurs  when  the  source- 
receiver  midpoint  is  at  the  eddy  center.  From  Eq.  (41a)  , it 
follows  that 


A4> 


N 


u>  B.U  r (R2  + B2N2  + B2R4N  2 + 2R2B^BJ1/2 
Lloo  2 3 23 

x {jQ[aR/(2ro) ] - 2JQ(a)}]  . 


(43) 


In  Eq.  (43),  there  are  no  current  effect^  as  all  currents  are 
perpendicular  to  the  transmission  range.  We  note  that  the 
position  of  the  source  and  the  receiver  are  not  unique,  since  they 
may  lie  on  any  eddy  diameter. 

We  first  choose  to  examine  the  range  dependence  of  the 
maximum  phase  shift  for  the  first  3 rays.  The  eddy  radius  and 


36 


1 


! 

:1 


1 


■ 


H 


maximum  current  speed  are  again  fixed  at  125  km  and  100  cm/sec ( 
respectively,  and  the  resulting  variations  with  range  illustrated 
in  Fig.  5.  The  phase  shift  for  each  ray  increases  from  R = 0,  where 
rays  are  reflected  d’ iectly  back  from  the  bottom.  In  this  case, 
the  AS’n  are  predicted  from  Eq.  (43)  to  be  linear  functions  of 
N,  so  the  phase  shifts  of  the  second  and  third  rays  are 
multiples  of  4.8  cycles,  the  shift  of  the  first.  We  observe 
from  Fig.  5 the  nonlinear  growth  of  the  phase  shifts  with 
range.  Moreover,  since  this  growth  changes  with  N,  there  is  a 
variation  in  the  relative  phase  spreading  among  the  rays, 
indicated  by  the  distances  between  curves  in  Fig.  5.  For  each 
positive  range,  the  A<t>N  are  nonlinear  functions  of  N so  that,  as 
range  increases,  the  relative  phase  spreading  between  the  second 
and  third  rays  becomes  significantly  less  than  that  of  the  first 
and  second.  For  example,  at  50  km,  the  phase  difference  between 
the  N=1  and  N=2  (N=2  and  N=3)  rays  is  4.0  (2.3)  cycles  compared 
to  3.0  (3.3)  cycles  at  R = 20  km. 

At  any  other  point  within  the  eddy,  similar  effects  would 
occur,  although  the  magnitudes  of  the  phase  variations  will  of 
course  be  smaller  than  the  maximum  values  appearing  in  Fig.  5. 

This  is  because  the  range-dependent  terms  in  Eq.  (41a)  are 
independent  of  source  and  receiver  position.  From  Eqs.  (41), 
the  current  effects  at  any  point  are  seen  to  be  linear  in  range  R,  so 
that  growth  of  current- i nduced  phase  effects  will  be  linear  in  R. 
Curves  illustrating  phase  variation  with  position,  such  as  those 
in  Fig.  3,  will  thus  grow  in  magnitude  with  increasing  range, 
although  the  relative  separation  will  vary  as  in  Fig.  5.  Any 
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asymmetries  resulting  from  current  effects  will  grow  linearly 
with  range  as  well. 

A striking  feature  of  a survey  of  eddy  observations,  such 
as  Ref.  1,  is  the  variability  encountered  in  size  and  strength. 

The  acoustical  consequences  of  such  variability  can  be 
investigated  by  changing  eddy  size,  described  by  the  radius  r , 
and  strength,  described  by  Uq.  We  first  fix  the  eddy  radius  at 
125  km  and  the  transmission  range  at  30  km,  and  study  the  effect 
of  variation  of  the  maximum  rotational  current  speed.  From 
Eq.  (43),  the  maximum  phase  variation  of  any  ray  is  a linear 
function  of  U . This  behavior  is  illustrated  in  Fig.  6 by 
both  the  magnitudes  of  per-ray  phase  shifts  and  relative  phese 
spreading  between  rays.  For  example,  at  Uq  = 150  cm/sec(  is 

20.2  cycles,  which  is  1.5  times  the  value  (13.5  cycles)  at 

UQ  = 100  cm/sec;  similarly,  A42  increases  by  a factor  of  1.5  to 

25.3  cycles.  At  other  source-receiver  locations  within  the  eddy, 
similar  effects  would  occur.  From  Eqs.  (41),  the  per-ray  phase 
shifts  resulting  from  both  currents  and  sound  speeds  are  linear 
in  Uq.  Hence,  curves  illustrating  phase  variations  during  an 
eddy  traversal,  such  as  those  in  Fig.  3,  would  increase 
similarly  as  functions  of  UQ,  in  both  magnitude  and  separation. 

The  effect  of  eddies  of  different  radii  can  be  investigated 
similarly,  by  fixing  Uq  = 100  cm/sec,  R = 30  km,  and  permitting 
rQ  to  vary.  The  A^  of  the  first  three  rays  for  radii 
characteristic  of  mesoscale  eddies,  illustrated  in  Fig.  7,  is 
seen  to  change  significantly  with  r . For  example,  for  the 
N=2  ray,  A4|2  changes  from  9.4  to  17.2  cycles  as  rQ  increases 
from  75  to  150  km.  From  Eq . (43),  A4>N  depends  on  rQ  both  ’''nearly 
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(as  with  U ) and  nonlinearly  throuqh  the  Bessel  function  J . 

o o 

However,  as  shown  in  Fig.  7 for  the  first  three  rays,  this 
nonlinearity  is  slight,  so  that  the  approximation 

A$n  = to  {B1UorQ(R2  + B2N2  + B2R4N2  + 2R2B2B3>1/2 

x [2  - 2 JQ(a) ] } (44) 

is  valid  for  small  aR/(2rQ).  We  note  that  this  expression  depends 
on  quantities  such  as  R and  UQ  in  the  same  way  as  Eqs.  (41a)  and 
(43),  so  it  could  be  used  to  generate  graphs  such  as  those  in 
Figs.  4 and  5. 

For  other  source-receiver  locations  within  the  eddy,  the 
term  in  Eq.  (41a)  due  to  sound-speed  effects  would  also  grow 
dramatically  with  increasing  eddy  radius.  In  contrast,  the  current- 
effect  term  does  not  have  a corresponding  factor  of  r . However, 
both  current  and  sound-speed  effects  on  per-ray  phase  depend  on 
the  relative  radii  at  the  source  (r  /r  ) and  receiver  (r  /r  ) 

o O HO 

in  Eqs.  (41).  Thus,  these  effects  at  a fixed  position  will  not 
be  linear  functions  of  r . For  example,  a source  that  is 
75  km  away  from  the  center  of  an  eddy  with  radius  of  75  km  is 
at  the  eddy  edge  (r  /r  = 1) , where  both  current  and  sound- 
speed  effects  are  absent.  On  the  other  hand,  in  an  eddy  of 
radius  150  km,  the  source  will  be  half-way  between  eddy  edge  and 
center  (r  /r  = 0.5),  where  currents  are  near  their  maxima,  and 

o O 

significant  sound-speed  effects  are  present.  Thus,  curves  such 
as  those  in  Fig.  3,  illustrating  phase  variations  during  an  eddy 
traversal,  would  vary  not  only  in  magnitude  and  relative 
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separation  but  also  in  horizontal  scale  as  eddy  radius  varies. 

In  summary,  Figs.  5,6,  and  7 illustrate  the  generality  of 

our  results,  as  a consequence  of  our  analytic  treatment  of  eddy 

environmental  and  acoustical  effects. 

V.  SUMMARY 

This  paper  represents  a general  investigation  of  the  effects 
of  mesoscale  cyclonic  eddies  of  arbitrary  size  and  strength  on 
relatively  short-range  underwater  sound  transmission.  Acoustically 
relevant  environmental  effects  are  determined  as  functions  of 
eddy  radius  and  maximum  current  speed,  through  use  of  a parametric 
eddy  model  which  specifies  fundamental  eddy  characteristics  (such 
as  SOFAR-axis  elevation  and  surface  currents  at  any  horizontal 
location  within  an  eddy).  Pivoting  bilinear  profiles  are  used 
to  describe  the  sound-speed  effects  of  an  eddy,  and  a linear 
approximation  to  the  current  is  made. 

Ray  geometry  for  a surfaced  source  and  receiver  within  an  eddy 
is  then  examined.  It  is  shown  that  the  three-dimensional  excursion 
from  the  vertical  source-receiver  plane  is  sufficiently  small 
that  rays  may  be  considered  as  two-dimensional.  Environmental 
discretization  of  eddy  sound-speed  and  current  effects  permits 
rays  to  be  approximated  by  a sequence  of  connected  circular  arcs, 
each  with  a modified  radius  of  curvature.  It  is  shown  that,  in 
general,  eddy  currents  can  cause  surface-ref lected/bot tom- 
reflected  rays  to  miss  a surface  or  bottom  reflection  and  become 
"mixed".  However,  none  of  those  rays  will  intersect  the 
receiver  for  ranges  of  less  than  50  km.  An  accurate  estimate 


for  the  launch  angle  of  all  rays  intersecting  the  receiver 
is  then  obtained,  with  effects  of  both  currents  and  horizontal 
sound-speed  variations  included.  The  geometry  of  all  rays 
is  then  completely  determined. 

The  travel  time  per  arrival  is  expressed  as  the  static 
travel  time  plus  small  changes  resulting  from  eddy  effects. 

These  changes  consist  of  terms  arising  from  both  sound  speed 
and  currents.  For  the  cold-core  cyclonic  eddies  considered  here, 
sound-speed  variations  result  in  a phase  delay  that  varies  with 
each  ray,  so  that  a relative  phase  spreading  of  arrivals  is 
observed.  The  effects  of  eddy  currents  can  result  in  either  a 
phase  advance  or  delay,  depending  on  the  current  relative  to  the 
source-receiver  direction.  The  current  effects  are  virtually 
independent  of  the  particular  rays,  so  that  a uniform  shift  in 
phase  is  observed. 

Several  examples  are  considered  to  illustrate  applications 
of  the  general  theory.  A fixed  source  and  receiver  separated 
by  30  km  are  moved  through  the  eddy,  and  the  resulting  phase 
shifts  of  the  principal  arrivals  computed.  It  is  observed  that 
both  the  individual  phase  shifts  and  the  relative  phase  spreading 
of  the  rays  reaches  a maximum  near  the  eddy  center,  where  sound- 
speed  perturbations  are  largest.  Current  effects  alone  are  shown 
to  produce  variations  of  about  one  cycle  between  identical 
positions  above  and  below  the  eddy  center.  This  is  a significant 
effect  when  compared  to  variations  induced  by  sound  speed.  The 
effect  of  source-receiver  orientation  within  the  eddy  is  considered 
next.  The  receiver  is  fixed  within  the  eddy,  and  the  source 
rotated  in  a circular  path  about  the  receiver.  Neglect  of  current 
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effects  is  shown  to  lead  to  an  error  as  large  as  20%  of  the  total 
variation  in  observed  phase.  Both  of  these  examples  indicate 
that,  in  eddy-acoustics  problems,  it  is  important  to  consider 
current  effects,  a factor  which  has  been  neglected  heretofore. 

General  results  are  presented  for  maximum  phase  fluctuation 
as  a function  of  transmission  range.  It  is  shown  that  changes 
in  this  quantity  grow  in  a nonlinear  manner  as  range  increases, 
but  that  relative  phase  spreading  is  not  a monotone  function  of 
range.  The  maximum  phase  variation  is  also  examined  as  a function 
of  eddy  size  and  strength,  to  illustrate  the  applicability  of 
our  general  procedure  to  arbitrary  eddies.  It  is  shown  that  phase 
fluctuations  of  over  20  cycles  can  occur  in  propagation  over  a 
30  km  range  within  an  eddy.  These  results  represent  an  initial 
step  by  the  authors  in  developing  general  results  applicable  to 
arbitrary  eddies  and  to  arbitrary  source  and  receiver  locations. 
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FIGURE 
FIG.  1 


FIG.  2 

FIG.  3 

FIG.  4 


FIG.  5 

FIG.  6. 

FIG.  7 


LEGENDS 

(a)  Eddy  sound-speed  effects  (solid  curves)  and 

bilinear  approximations  (dashed  lines);  Uq  = 100  cm/sec, 

r = 125  km.  (b)  An  eddy  current  (solid  curve)  and 
o 

linear  approximation  (dashed  line) . 

(a)  Source  and  receiver  coordinates.  (b)  Idealized 

ray,  showing  reflection  points  x^  and  environmental 

discretization  intervals  I.. 

i 

Phase  shifts  of  first  three  SRBR  arrivals  when  trans- 
mission range  varies  as  shown;  UQ  = 100  cm/sec, 
rQ  = 125  km,  R = 30  km,  f = 400  Hz. 

Phase  shift  of  N = 1 ray  when  transmission  range  rotates 
through  the  angle  8 as  shown.  Solid  curve  indicates 
currents  present;  dashed  curve  neglects  currents; 

UQ  = 100  cm/sec,  rQ  = 125  km,  R = 30  km,  f = 400  Hz. 

Maximum  change  in  phase  of  first  three  arrivals  as  a 
function  of  transmission  range;  UQ  = 100  cm/sec, 
rQ  = 125  km,  f = 400  Hz. 

Maximum  change  in  phase  of  first  three  arrivals  as  a 
function  of  eddy  strength;  rQ  = 125  km,  R = 30  km, 
f = 400  Hz. 

Maximum  change  in  phase  of  first  three  arrivals  as  a 
function  of  eddy  radius;  UQ  = 100  cm/sec,  R = 30  km, 
f = 400  Hz. 
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